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The PHENIX experiment at RHIC has measured charged hadron yields at mid-rapidity over a 
wide range of transverse momentum (0.5 < pr < 10 GeV/c) in Au + Au collisions at ^/sNN = 200 
GeV. The data are compared to tt" measurements from the same experiment. For both charged 
hadrons and neutral pions, the yields per nucleon-nucleon collision are significantly suppressed in 
central compared to peripheral and nucleon-nucleon collisions. The suppression sets in gradually 
and increases with increasing centrality of the collisions. Above 4-5 GeV/c in pr, a constant and 
almost identical suppression of charged hadrons and 7r'''s is observed. The px spectra are compared 
to published spectra from Au -\- Au at y^ijvjv = 130 in terms of xt scaling. Central and peripheral 
7r° as well as peripheral charged spectra exhibit the same xt scaling as observed in p -I- p data. 

PACS numbers: 25.75.Dw 
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I. INTRODUCTION 

Lattice Quantum-Chromo-Dynamic (QCD) calcula- 
tions predict a new state of matter of deconfined 
quarks and gluons at an energy density exceeding 
~1 GeV/fm^ It has long been suggested that such 
a "quark gluon plasma" may be produced in collisions 
between ultra-relativistic heavy nuclei Indeed, mea- 
surements of transverse energy produced in high energy 
Pb + Pb and Au + Au collision suggest that energy den- 
sities above 3 GeV/fm^ at SPS pf and 5 GeV/fm^ at 
RHIC 0, 0] have been reached. However, this conclu- 
sion relies on model assumptions |E S IS 13 to relate the 
properties of the hadronic final state to the initial state 
dynamics. 

The spectra of high transverse momentum {pt) 
hadrons resulting from the fragmentation of hard- 
scattered partons potentially provide a direct probe of 
the properties of the initial state. Theoretical calcula- 
tions show that the outgoing high-pT partons radiate sub- 
stantially more energy when propagating through dense 
matter than when propagating in the vacuum, resulting 
in a softening of the hadron px spectrum [iDj , with the 
energy loss of the partons depending on the gluon den- 
sity of the matter [lllll^ . Formation time considerations 
suggest that hard scattered partons are "produced" at 
the earliest stage of the collision, thus directly probe the 
dense matter from the time of their creation. Therefore, 
a detailed analysis of high-pT hadron production may re- 
veal information on the properties of the dense medium 
created early in the collisions 0, 0, 0| . 

At the energies reached at the Relativistic Heavy Ion 
Collider (RHIC), high-p^ hadrons are copiously pro- 
duced. In nucleon-nucleon collisions, it has been well 
established that hadrons with pr > "2 GeV/c result pri- 
marily from the fragmentation of hard-scattered partons, 
and that the px spectra of these hadrons can be calcu- 
lated using perturbative QCD (pQCD) Initial 
measurements of hadron pT spectra in Au + Au collisions 
at y^SNN = 130 GeV led to the discovery of a substantial 
suppression of hadron yields per nucleon-nucleon colli- 
sion relative to pp data 0, 0, Holl. Data from ^/snn 
= 200 GeV confirm these resultsU^, HI IH E^. The 
suppression is observed in central but not in peripheral 
collisions. These observations are consistent with pQCD- 
inspired modelling of parton energy loss in dense mat- 
ter 0il23- However, alternative interpretations that do 
not assume the formation of a deconfined phase have been 
proposed based on the modifications of the parton dis- 
tribution functions in the initial state or final-state 
hadronic interactions [27| . 

In addition to hadron suppression, an unexpectedly 
large fraction of baryons has been observed in central 
Au + Au collisions for pr up to 4-5 GeV/c pll29ll30j. 
which complicates the interpretation of the high px 
results. The observed baryon to meson ratio from 



PHENIX is inconsistent with jet fragmentation in 
p + p and e+e~ colhsions [s^l- While the origin of 
this effect is unclear, it could point towards bulk particle 
production ( "soft physics" ) contributing to the px spec- 
tra out to 4-5 GeV/c. It has been suggested that coales- 
cence of thermalized quarks combining with energy loss 
of hard-scattered partons can account for the unusual 
particle composition, which shifts the region dominated 
by hard-scattering to higher px [13 ■ 

Systematic measurements of the px, centrality, parti- 
cle species, and ^'snn dependence of the suppression can 
constrain competing descriptions of high-pr hadron pro- 
duction. In this paper, we present new data on inclusive 
charged hadron production for 0.5 < pt < 10 GeV/c, 
measured over a broad range of centrality in Au + Au 
collisions at ^s^vw = 200 GeV by the PHENIX Collab- 
oration at RHIC. These data are compared to data on 
neutral pion production |2^ and to data from Au + Au 
collisions at ^/snn = 130 GeV 0, JJ^, all measured 
within the same experiment. 

The remainder of the paper is organized as follows. 
Section Ogives a detailed account of the charged parti- 
cle analysis. Centrality and px dependence of the charged 
hadron px spectra are discussed in Section IIII Al Sec- 
tion 1111 ij| studies the charged hadron suppression and 
compares the results to tt" data. In Section IIII CI we dis- 
cuss the y/SNN dependence of both charged hadron and 
neutral pion production and test possible a;T-scaling. A 
summary is given in Section llVI 



II. DATA ANALYSIS 
A. PHENIX Detector 

The PHENIX experiment consists of four spectrome- 
ter arms — two around mid-rapidity (the central arms) 
and two at forward rapidity (the muon arms) — and a 
set of global detectors. The central arm and south Muon 
arm detectors were completed in 2001 and took data dur- 
ing Au + Au operation of RHIC the same year (RUN-2). 
The layout of the PHENIX experiment during RUN-2 is 
shown in Figure^ Each central arm covers the pseudo- 
rapidity range \r]\ < 0.35 and 90 degrees in azimuthal 
angle (p. In each of the central arms, charged particles 
are tracked by a drift chamber (DC) positioned from 2.0 
to 2.4m radially outward from the beam axis and 2 or 
3 layers of pixel pad chambers (PCI, (PC2), PC3 lo- 
cated at 2.4m, (4.2m), 5m in radial direction, respec- 
tively). Particle identification is provided by ring imag- 
ing Cerenkov counters (RICH), a time of flight scintillator 
wall (TOF), and two types of electromagnetic calorime- 
ters (lead scintillator (PBSC) and lead glass (PBGL)). 
The magnetic field for the central spectrometer is axially 
symmetric around the beam axis. Its component par- 
allel to the beam axis has an approximately Gaussian 
dependence on the radial distance from the beam axis, 
dropping from 0.48 T at the center to 0.096 T (0.048 T) 
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FIG. 1: (Color online) PHENIX experimental layout for the 
Au + Au run in 2001. The top panel shows the PHENIX 
central arm spectrometers viewed along the beam axis. The 
bottom panel shows a side view of the PHENIX muon arm 
spectrometers. 

at the inner (outer) radius of the DC. A pair of Zero- 
Degree Calorimeters (ZDC) and a pair of Beam-Beam 
Counters (BBC) were used for global event characteriza- 
tion. Further details about the design and performance 
of PHENIX can be found in |3|. 



B. Event Selection 

During RUN-2, PHENIX sampled an integrated lumi- 
nosity of 24 fj,b~^ for Au + Au collisions at ^/snn — 200 
GeV. Minimum bias events were selected by a coincidence 
between the ZDCs and the BBCs. This selection corre- 
sponds to 92.2^3 0% of the 6.9 barn Au + Au inelastic 
cross section. The event centrality is determined by cor- 
relating the charge detected in the BBCs with the energy 
measured in the ZDCs. Two sets of centrality definitions 
are used in this analysis: a "Fine" set of centralities, 
which corresponds to 0-5%, ...,15-20%,20-30%,.. .,80-92%, 
and a "Coarse^^ set of centralities, which corresponds 
to 0-10%,10-20%,20-30%^,80-92%. A Glauber model 
Monte-Carlo simulation [ssIIb^ . l3^ l3^ that includes the 
responses of BBC and ZDC gives an estimate of the av- 



TABLE L Centrality classes, average number of A'^ -I- A'' col- 
lisions, average number of participant nucleons, and average 
nuclear overlap function obtained from a Glauber Monte- 
Carlo simulation of the BBC and ZDC responses for Au + Au 
at y^sjvjv ~ 200 GeV. Each centrality class is expressed as a 
percentage of (tauAu = 6.9 b. Two sets of centrality defini- 
tions are used in this analysis: a "Fine" set of centralities, 
which corresponds to 0-5%, ...,15-20%,20-30%,... ,80-92%, and 
a "Coarse" set of centralities, which corresponds to 0-10%, 10- 
20%,20-30%,... ,80-92%. 



Centrality 




{Npart ) 


{TAuAu){mb~'') 


- 


5% 


1065 ± 105.5 


351.4 ±2.9 


25.37 ± 1.77 


5 - 


10% 


854.4 ±82.1 


299 ± 3.8 


20.13 ± 1.36 


10 - 


15% 


672.4 ± 66.8 


253.9 ±4.3 


16.01 ± 1.15 


15 - 


20% 


532.7 ± 52.1 


215.3 ±5.3 


12.68 ±0.86 


- 


10% 


955.4 ± 93.6 


325.2 ±3.3 


22.75 ± 1.56 


10 - 


20% 


602.6 ± 59.3 


234.6 ±4.7 


14.35 ± 1.00 


20 - 


30% 


373.8 ± 39.6 


166.6 ±5.4 


8.90 ± 0.72 


30 - 


40% 


219.8 ±22.6 


114.2 ±4.4 


5.23 ±0.44 


40 - 


50% 


120.3 ± 13.7 


74.4 ±3.8 


2.86 ±0.28 


50 - 


60% 


61.0 ±9.9 


45.5 ±3.3 


1.45 ±0.23 


60 - 


70% 


28.5 ± 7.6 


25.7 ±3.8 


0.68 ±0.18 


70 - 


80% 


12.4 ±4.2 


13.4 ±3.0 


0.30 ±0.10 


80 - 


92% 


4.9 ± 1.2 


6.3 ± 1.2 


0.12 ±0.03 


60 - 


92% 


14.5 ±4 


14.5 ±2.5 


0.35 ±0.10 


min 


bias 


257.8 ±25.4 


109.1 ±4.1 


6.14 ±0.45 



erage number of binary collisions (Ncou), participating 
nucleons (Npart) and nuclear overlap function (TauAu) 
for each centrality class. The calculated values of (Ncoii), 
(Npart) and (TauAu) for each centrality class are listed in 
Table n 

In addition to the event selection, the BBCs also allow 
us to reconstruct the collision vertex in the beam direc- 
tion (z) with a resolution of 0.5 cm. An offline z-vertex 
cut, \zytx\ <30 cm, was applied to the minimum bias 
events. After this selection, a total of 27 x 10® minimum 
bias Au + Au events were analyzed to obtain the charged 
hadron spectra presented in this paper. 

C. Charged Particle Tracking and Momentum 
Measurement 

Charged hadron tracks are measured using informa- 
tion from the DC, PCI, PC2 and PCS detectors of the 
west central-arm and the BBC. The projections of the 
charged particle trajectories into a plane perpendicular 
to the beam axis are detected typically in 12 wire planes 
in the DC. The wire planes spaced at 0.6 cm intervals 
along the radial direction from the beam axis. Each 
wire provides a projective measurement, with better than 
150 fim spacial resolution in the azimuthal {(/)) direction. 
Eight additional wire planes in the DC provide stereo- 
scopic projections, which together with the space point 
measured at the PCI and the vertex position measured 
by the BBC determine the polar angle of the track. Tra- 
jectories are confirmed by requiring matching hits at both 
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PC2 and PCS to reduce the secondary background. 

Tracks are then projected back to the coUision vertex 
through the magnetic field to determine the momentum 
p. The transverse momentum pT is related to the de- 
flection angle a measured at the DC with respect to an 
infinite momentum trajectory. For tracks emitted per- 
pendicular to the beam axis, this relation can be approx- 
imated by 



K 

PT 



(1) 



where if = 87 mrad GeV/c is the effective field integral. 

The momentum scale is verified by comparing the 
known proton mass to the value measured for charged 
particles identified as protons from their time-of-flight. 
The fiight-time is measured in the TOF detector, which 
cover 7r/4 of the azimuthal acceptance in the east arm. 
The absolute value of the momentum scale is known to 
be correct to better than 0.7%. 

The momentum resolution is directly related to the a 
resolution, 



5p/p ~ 5a I a 



(2) 



K 



where 5a is the measured angular spread, which can be 
decomposed into the contribution from multiple scatter- 
ing {(Jms) and the contribution from the intrinsic point- 
ing resolution (ctq) of the DC. At high py, ctq is the 
dominating contribution, i.e. 5a — ctq. We measure 
ctq, ~ 0.84 ± 0.05 mrad/(GeV/c) using zero field data, 
where we select high-momentum tracks by requiring en- 
ergetic hadronic showers in the electromagnetic calorime- 
ters. The width of the proton mass as function of p^ 
independently confirms the momentum resolution. In 
summary, the momentum resolution is determined to 
be Sp/p ~ 0.7% © 1.0% p (GeV/c). Further details on 
track reconstruction and momentum determination can 
be found in l33l. 



D. Background Rejection and Subtraction 

Approximately 95% of the tracks reconstructed by the 
DC originate from the event vertex. The remainder have 
to be investigated as potential background to the charged 
particle measurement. The main background sources in- 
clude secondary particles from decays and e^e~ pairs 
from the conversion of photons in materials between the 
vertex and the DC. Depending on how close the conver- 
sion or decay point is to the DC, or depending on the 
Q-value of the decay, these tracks may have a small de- 
flection angle a at the DC. Thus, according to Eq.^] they 
are incorrectly assigned a large momentum. In this anal- 
ysis, the pt range over which charged particle production 
is accessible in PHENIX is limited by this background. 
We exploit the track match to PC2 and PCS to reject as 



much of the background as possible, then employ a sta- 
tistical method to measure and subtract the irreducible 
background. 

For primary tracks, the distance in both the r-0 and 
the z direction between the track projection point and 
the measured PC hit position is approximately Gaussian 
with a mean of and a width given by. 



^raatch 



' det 



'itch 



pP 



(3) 



where cr^J;"*'^'' is the finite detector resolution (which in- 
cludes DC pointing (or a) resolution and the PC2, PCS 
spacial resolution) , and cr™"*'^'* is the multiple scattering 
contribution. 

Despite being incorrectly reconstructed with large px, 
the majority of the background particles have low mo- 
menta. While travelling from the DC to the PC2 and 
PCS, they multiple scatter and receive an additional de- 
flection from the fringe field. This causes a correlated 
defiection between the measured positions at PC2, PCS, 
and the projections calculated from tracks measured by 
the DC and PCI. The displacements in r-(/) and z direc- 
tions are represented by and D^. Since the residual 
bend depends on the z component of the fringe field, 
which decreases rapidly at large |ry|, a fiducial cut of 
I77I <0.18 was applied to ensure that the residual bend 
due to the fringe field is almost independent of z. We 
focus on the displacement in T-(f), D^, which are large 
for low momentum tracks due to the residual bend. The 
ZJ^'s at PC2 and PCS are correlated with each other, as 
shown in Figure|21 Most of the tracks lie in a narrow win- 
dow around the diagonal line. The width of this window 
is given by the PC2 and PCS detector resolutions, which 
are of the order of a few millimeters. Multiple scatter- 
ing and residual bend broaden the matching distribution 
along the diagonal line. To optimize background rejec- 
tion, we define two orthogonal projections. 



(4) 



D'^ is the variable along the correlated direction, is 
the direction normal to . A ±2tT cut on these variables 
is applied in the data analysis. In the remaining discus- 
sion, unless stated otherwise, only tracks satisfying these 
cuts are included. 

After matching cuts, the background level is less than 
6% for Pt < 4: GeV/c, but increases rapidly at higher pT- 
For 4 < Pt < 10 GeV/c, the most significant remaining 
background sources are e~^e~ from conversion of pho- 
tons close to the DC and particles from weak decays of 
long lived particles, mainly of and if ° . These back- 
grounds are estimated and subtracted separately from 
the D'^ distribution for all tracks, as will be discussed in 
the rest of this section. 
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FIG. 2: -D^*^^ (the difference between projection and iiit lo- 
cation in r-ff) direction at PC2) versus D^'^'' in centimeters for 
tracks with reconstructed pt > 4: GeV/c. PC2, PCS matching 
differences are correlated, with signal tracks peaked around 
and background tracks extend along the direction. The 

double-peak structure along is related to the finite gran- 
ularity of PC2 and PCS pads. The positive directions of D'^ 
and are indicated by the arrow. A ±2a cut on these 
variables is illustrated by the box region inside the dashed 
lines. 



To separate the two background sources, we take ad- 
vantage of the RICH to tag electrons. Charged parti- 
cles with velocities above the Cherenkov threshold jth = 
35 (CO2 radiator) will emit Cherenkov photons, which 
are detected by photon multiplier tubes (PMT) in the 
RICH We characterize the Cherenkov photon yield 
for a charged particle by Npmt, the number of PMTs 
with signals above threshold associated to the track. For 
reconstructed electrons {pr > 150MeV/c), the average 
number of associated PMTs is (Npmt) ~ 4.5. The prob- 
ability to find at least one PMT above threshold is more 
than 99%. For pions, the Cherenkov threshold is 4.8 
GeV/c, and the number of associated PMTs reaches its 
asymptotic value only well above 10 GeV/c; {Npmt) in- 
creases from 1.4 at 6 GeV/c to 2.8 at 8 GeV/c and 3.6 
at 10 GeV/c. 

Tracks (Np) with at least one associated RICH PMT 
contain both conversion electrons and real pions. Their 
matching distributions in D'^ are presented in Figure |3| 
for a sample range of 6 < pt < 7 GeV/c. Also shown is 
the matching distribution for conversion electrons from 
Monte-Carlo simulation. The contributions from pions 
and electrons are clearly distinguishable. For pions with 
Pt < 10 GeV/c, {Npmt} has not reached its asymptotic 
value. A requirement of Npmt > 5 rejects most of the 
pions while preserving a well-defined fraction (Re) of the 
electrons. To measure Re from the data, we select tracks 



with an apparent px > 10 GeV/c. ^ The fraction of 
tracks with Npmt > 5 is measured to be Re — 0.458 ± 
0.05. Both Monte-Carlo and data show a small variation 
of Re with Pt and centrality. This variation is included 
in the error on Re- The total electron background is 
calculated using tracks with Npmt > 5 (iVg) as, Ne/Re- 
The number of real pions in the RICH-associated sample 
for each pT bin is then calculated as, 

Sr = Nr-^ . (5) 

Ke 

With this method, a small fraction of genuine pions, 
which satisfy Npmt > 5, is subtracted. This fraction is 
negligible below 7 GeV/c, but increases rapidly towards 
higher px- This loss is corrected using the PHENIX 
Monte-Carlo simulation. In this 100% error on 

the correction is assigned. 

The sample of tracks {N^p) with no associated RICH 
PMT contains a mixture of tt, K,p, contaminated by the 
decay background. Their matching distributions in 
are shown in Figure 01 for 6 < pt < 7 GeV/c, together 
with the matching distribution for decay particles from 
MC simulation. A Monte-Carlo study shows that the ap- 
parent momentum of these tracks is nearly uncorrelated 
with true momentum and therefore the distribution of 
this background in is nearly independent of the ap- 
parent momentum. We select a nearly pure background 
sample using tracks with reconstructed px > 10 GeV/c 
and measure the ratio of the number of tracks passing a 
\D'^\ < 2(7 cut to tracks in the interval 3a < \D'^ \ < 9a: 

^ ^ Nr,R{pT>lOGeV/c,\D+\<2a) 

Nnr{pt > IQGeV/c, 3a < \D+\ < 9a) ^ ' 
= 0.424 ±0.05 

The error quoted takes into account the small variation 
of Rdecay with px and centrality. Since the average yield 
of real hadrons in this interval is small, we estimate the 
decay contribution as a function of px to be NNp{3a < 
< 9a) X Rdecay Finally, the signal is calculated as, 

Snr = NNR{\D+\<2a)^ (7) 
NNR{3a <\D+\< 9a) x Rd ecay 

Figure gives the total signal, obtained as Sr + Snr, 
with the decay and conversion background subtracted. 
On the right hand side, the signal-to-background ratio is 
shown. The background increases with increasing pT • At 



In this px range, the background yield decreases slowly with 
px, while the true tt yield decreases rapidly as px incr ease s. By 
comparing the measured tt'' spectrum from PHENIX |21| with 
the charged hadron spectrum before background subtraction at 
Px > 10 GeV/c, the integrated signal yield is estimated to be 
less than 3% and thus may be neglected. 
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FIG. 3: (Color online) Background contamination due to electrons, illustrated by the track match in for tracks with 
associated RICH PMTs and 6 < pr < 7 GeV/c. The matching distributions are shown for minimum bias events and separately 
for positive (left) and negative (right) charged tracks. The first three distributions represent the raw counts for all tracks 
with RICH association (thick solid line), estimated conversion backgrounds (dashed line) and charged pions (dot-dashed line) 
that were obtained by subtracting the dashed line from the solid line. The thin solid line represents the matching distribution 
of background electrons from Monte-carlo simulation, arbitrarily scaled to match the data. The 2 a matching window are 
illustrated by the vertical dashed line. Since e+ and e~ are deflected in opposite directions by the fringe field, they are shifted 
to positive and negative directions, respectively. 




FIG. 4: (Color online) Background contamination due to decays, illustrated by the track match in for tracks without an 
associated RICH PMT and with 6 < pr < 7 GeV/c, shown for minimum bias events and separately for positive (left) and 
negative (right) charged tracks. The flrst three distributions represent the raw counts for all tracks without a RICH association 
(thick solid line), estimated decay backgrounds (dashed line), and signal tracks (dot-dashed line) that were calculated as the 
difference of the two. The thin solid line represents the matching distribution of decay background from Monte-carlo simulation, 
arbitrarily scaled to match the data. The 2 a matching window are illustrated by the vertical dashed line. Outside the signal 
window, the shape of the dashed line matches the solid line rather well, the difference of 10% level is taken into account in the 
error estimation of Rdecay 



4 GeV/c the signal-to-background ratio is about 10, and 
decreases to 1 at 7.5 GeV/c and to ^0.3 at 10 GeV/c. 

Weak decays of short hved particles, mainly K^, A 
and A within the magnetic field provide an additional 
source of background. A significant fraction of this back- 
ground is subtracted using the Rdecay method described 
above. However, secondary particles from decays close 
to the event vertex are not subtracted since they are 
nearly indistinguishable from primary particles. This 
"feed-down" contaminates the track sample without the 



associated RICH PMTs, Snr (about 40% of all charged 
particles at high pt), and needs to be subtracted from 
the data. 

To estimate the feed-down contribution we generate 
Au + Au events with HIJING 41J , reconstruct them 
through the PHENIX Monte-Carlo simulation, and count 
the secondaries which survive all analysis cuts. The sec- 
ondaries from A and A decays are counted relative to the 
reconstructed (p-f p) , and correspondingly, those from 
relative to {K+ + if")/2. We tune the (A -f A)/(p -t- p). 
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FIG. 5: (Color online) Amount of background estimated as function of px for minimum bias collisions. The left figure shows 
the background subtracted charged hadron spectra (filled square) , the background from 7 conversions (open square) and decays 
(open triangle). The right figure shows the signal to background ratio. Only statistical errors are shown. 



K'^ / {0.5{K^ + K^)) ratios by weighting the particle dis- 
tributions generated according to HIJING such that they 
reproduce the nearly py-independent experimentally ob- 
served ratios from Au + Au colHsions at ^/snn — 130 
GeV lillil. 



The final feed-down contribution depends on the choice 
of the A and pT spectra and of their yields in the 
high PT range where they are not measured. Both yields 
and spectral shapes are varied within limits imposed by 
the spectrum for tracks that do not fire the RICH. The 
average feed-down contribution depends on px and varies 
between 6 to 13% relative to the total charged hadron 
yield; it is subtracted from the charged spectra. The 
systematic uncertainties are estimated from the spread 
of the feed-down contributions obtained with different 
assumptions. The uncertainties are approximately 60% 
of the subtraction, and depend on pT and centrality. 

Tablenisummarizes the systematic errors on the back- 
ground subtraction ^. All errors are correlated with px 
and are presented as relative errors to the charged hadron 
yield. The uncertainty of the pion oversubtraction correc- 
tion (SttIoss) was re-scaled by the fraction of signal tracks 
with RICH association, i.e. Su/{Su + Snr)- Errors on 
the scaling factors Re and Rdecay were individually folded 
with the signal-to-background ratios in the two samples. 
The resulting uncertainties on the charged yields were 
then added in quadrature (<5fl^©_R,jecoH)- uncertainty 
of the K'^, a, and A feed-down subtraction is denoted by 

^ f eeddown ■ 



^ We should emphasize that, in the remaining discussion unless 
stated otherwise, all systematic errors listed in Tables have been 
adjusted to Icr errors. 



E. Corrections and systematic uncertainties 

After background subtraction, we have determined a 
single, pT-dependent correction function to correct the 
hadron spectra for acceptance, decay in flight, recon- 
struction efficiency and momentum resolution. This 
function is determined using a GEANT Monte-Carlo 
simulation 0| of the PHENIX detector in which sim- 
ulated single tracks are reconstructed using the same 
analysis chain applied to the real data. Because of de- 
cays and multiple scattering, the correction function de- 
pends on the particle species. This is reflected in Fig- 
ure El where the correction functions averaged between 
TT"*" and 7r~, and K~ , p^ and p~ are shown sepa- 
rately. At pt < 3 GeV/c, the kaon correction function is 
significantly larger than those for pions and protons; For 
Pt > 3 GeV/c, this difference is less than 15%. To take 
into account this species dependence, we determine the 
correction function separately for pions, kaons, protons, 
and their anti-particles. The final correction function is 
then obtained by combining the correction functions for 
the different particle species weighted by the measured 
PT-dependent particle composition from |4^ . Above 2 
GeV/c, where kaon data are not available, we assume 
the K/tt ratio is constant within ±10% from the value 
observed at 2 GeV/c. This assumption leads to a 2.5% 
systematic error in the correction function. The resulting 
correction function is plotted in the upper left panel of 
Figure [3 The sharp rise below 2 GeV/c is due to loss 
in acceptance and decays in flight. Above 2 GeV/c, the 
correction decreases only slowly with pT- For > 4 
GeV/c, the correction varies by less than ±5%. 

The data are also corrected for efficiency losses due to 
detector occupancy. Though this is negligible for periph- 
eral collisions, these losses are important in central colli- 
sions, and are evaluated by embedding simulated tracks 
into real events. The average track reconstruction ef- 
ficiency in the active detector area is larger than 98% 
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TABLE II: Systematic errors on background subtraction. All errors are given in percent and are quoted as la errors. These 
errors are correlated with pt- 



Pt (GeV/c) S^ioss (%) Sii^(DRdacayi%) Sf^eddown (%) total(%) 

< 5 0.3 0.3 5 5 

5 - 6 0.6 1.8 5 5.3 

6 - 7 1.4 4.1 8.5 9.5 

7 - 8 4.6 7.1 7.8 11.5 

8 - 9 9.9 17.6 6 21.1 

9 - 10 19.4 23.5 6 31.1 
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FIG. 6: (Color online) Averaged correction functions for tt'*' 
and TT^, p and p, and and K~ . 



for peripheral collisions, but decreases to 70 ± 3.5% for 
central collisions. As shown in the lower part of Fig- 
ure [71 the efficiency loss is independent of px within 
a ±3% systematic uncertainty from 1.5 to 10 GcV/c. 
Based on this observation, the full correction can be 
factorized into centrality-dependent (i.e. detector oc- 
cupancy dependent) correction function, c{Npart), and 
PT-dependent correction function, c{pt)- The centrality- 
dependent correction function is shown on the upper 
right panel of Figure Most of the efficiency loss is 
due to hit overlaps, which can shift the hit positions in 
the DC or PC's outside of the matching windows. The 
±2(T matching windows are larger at low pT to account 
for multiple scattering (see Eq.|3Jl, thus the tracks are less 
vulnerable to the effect of hit merging. This effect has 
been taken into account by applying a slightly smaller, 
PT-dependent, occupancy correction at px < 1-5 GeV/c. 

Figurc0also shows systematic errors on the correction 
functions. These errors include not only the errors on 
the correction itself, but also the uncertainty due to the 
background subtraction procedure. 

Finally, the inclusive charged hadron yield are ob- 
tained by multiplying the pT-dependent correction func- 
tion, c{pt), and centrality-dependent correction function. 



c(Npart), with the background subtracted spectra and di- 
viding by the number of events for every centrality class 
as: 



1 



dN 



1 



Ne^ts 2TTpTdpTdri 



X c{pt) X c{NpaTt) (8) 



dN 



PrdpTdrj 



hgr ~ subtracted 



The systematic errors on the spectra, which are com- 
mon to all centrality classes, are listed in Table Hn 
Sources of systematic uncertainties are: the matching 
cuts {5match), normalization (<5„orm), particle composi- 
tion {5mix)i momentum resolution (<5reso)i momentum 
scale {5 scale) ^ and background subtraction (Sbgr) from Ta- 
ble^ The normalization error is independent of Pt- All 
other errors vary with px but are highly correlated bin- 
to-bin, which means that points in neighboring pT bins 
can move in the same direction by similar factors. 

The centrality-dependent systematic errors are quan- 
tified in terms of the central-to-peripheral ratio. Rep, as 
given in Table IIVI Besides the uncertainty on the oc- 
cupancy correction {S occupancy) illustrated in Figure [7| 
the background subtraction procedure has a centrality- 
dependent uncertainty. As discussed in Section III Dl the 
errors on Re and Rdecay reflect the pr and centrality de- 
pendencies. The centrality-dependent part contributes 
about half of the error on both i?e and Rdecay , and hence 
does not cancel in Rep. Since the errors on R^ and Rdecay 
are independent, the uncertainty on R^p is approximately 
equal to 5r^(dr^^^^^ from Table HIl Finally, Sfeeddown is 
the centrality-dependent error from feed-down subtrac- 
tion. 



III. RESULTS 

A. Inclusive charged hadron pr spectra 

Figure|Slshows the inclusive charged hadron pt spectra 
for various centrality classes. All spectra exhibit power- 
law tails at high pT- But for peripheral collisions, the 
power-law shape is more concave than for central colli- 
sions. More details of the centrality dependence of the 
spectral shape can be seen from Figure El which shows 
for each centrality class the ratio of the spectra to the 



10 



c 

u 
c 

3 

H— 

c 
o 

u 

0) 

1_ 
1_ 

o 
o 





3456789 10 

PT(GeV/c) 



50 100 150 200 250 300 350 400 

Npar, 



g 

(0 




8 9 10 

Pp(GeV/c) 

FIG. 7: (Color online) Functions used to correct the charged particle pr spectra. The upper left panel shows the pr dependent 
correction, c{pt)- The upper right panel shows the centrality dependent correction, c{Npart)- Systematic uncertainties are 
indicated by the dashed lines. The two corrections factorize at pr > 1.5 GeV/c, so that for given centrality the full correction 
function is given by c(pt) x c{Npart)- The accuracy of this factorization is demonstrated in the lower panel. The ratio of the 
full correction for central collisions (5% most central) to the correction for single particle events varies by less than 3% above 
1.5 GeV/c (the error bar is the statistical error from the Monte-Carlo calculation). 



TABLE III: 


Systematic errors 


on the pt spectra. 


All errors are 


given in percent and 


are quoted as 


1(7 errors. 


They are either 


normalization errors or are pr 


correlated errors. 












PT (GeV/c) 


Smatch{%) 


^norm (9o) 






Sscale{%) 
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< 1 
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0.6 


0.6 
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0.6 
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0.6 


3.1 
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8 - 9 
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3.2 


1.8 
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3.1 


21.1 


21.9 


9 - 10 
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3.2 


1.8 


5.3 
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minimum-bias spectrum. In these ratios, most system- 
atic errors cancel or affect the overall scale only. The 
characteristic centrality dependence of the shape alre ady 
observed in ^snn = 130 GeV Au + Au collisions \\A\\m 
is more apparent at ^Jsnn = 200 GeV. In peripheral col- 
lisions, the ratio decreases up to ^ 2 GeV/c and then 
rises up to about 4 GeV/c. The trends are reversed in the 
most central collisions. In the range above 4-5 GeV/c, 
all ratios appear to be constant as function of pT^ which 
would imply that they have a similar centrality indepen- 
dent shape. 

Based on the different trends observed in Figure |51 
we can distinguish three pt regions: 0.5-2, 2-4.5 and 
> 4.5 GeV/c. The different centrality dependence of the 



spectral shape in these regions can be quantified by a 
truncated average pr'- 



8 GeV/c 



Pt • dN/ dpT 



jp. 



8 GeV/c 



dN/dpi 



mm 

' Pt I 



(9) 



which is insensitive to the normalization of the spectra. 
The upper bound of 8 GeV/c in the integral is given by 
the limited pt reach for peripheral centrality classes as 
shown in Figure|H| In Figure^! the values of (p^""'^) for 
the three p™™ values are plotted as function of central- 
ity, represented by the average number of participating 
nucleons {Npart) for each centrality class. 
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TABLE IV: Systematic errors on the central-to-peripheral ratio. All errors are given in percent and are quoted as la errors. 
Most of the systematic errors listed in Table IITTI cancel in the central-to-peripheral ratio. Only those errors that are uncorrelated 
with centrality are shown here. 



PT (GeV/c) 5occupancy{%) 5feeddon,n{%) SB.^(SRdacav (-^) total(%) 

< 6 5 5 1.8 7.3 

6 - 7 5 5 4.1 8.2 

7 - 8 5 5 7.1 10 

8 - 9 5 5 17.6 19 

9 - 10 5 5 23.5 24.6 
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FIG. 8: (Color online) pr spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes 
derived from the pseudo-rapidity region \rj\ < 0.18. The minimum bias spectrum has been multiplied by 5 for visibility. Only 
statistical errors are shown in the spectra. Most of the pr dependent systematic errors are independent of centrality and are 
tabulated in Table ITTTl 



For 

pmin _ Q g GeV/c, where particle production is ex- 
pected to be governed by soft physics, (p^""'^) increases 
with Npart- This trend is also seen for the average pt 
of identified charged hadrons, and reflects the increased 
radial flow of soft particles in more central collisions . 
For p™™ = 2 GeV/c, the trend is significantly differ- 
ent. For peripheral collisions, (p^"""^) is substantially 
larger than the value obtained with p™*" = 0.5 GeV/c 
due to the presence of the power-law tail. With increas- 
ing Npart, for p^" = 2 GeV/c decreases and the 
values obtained with p™™ — 0.5 and 2 GeV/c approach 
each other, which indicates an almost exponential spec- 
trum in central collisions between 0.5 and 2 GeV/c. For 
the highest px range (p™™ = 4.5 GeV/c), {ptf^"'^) is ap- 
proximately constant. This implies that the shape of the 
spectrum is nearly independent of centrality, as would be 
expected if this region is dominated by hard scattering. 

However, the yields at high pT do not scale with the 
number of nucleon-nucleon collisions; they are suppressed 
comparing to the binary collision scaling expected for 
hard scattering processes. This can be clearly seen from 
Figure [TTl which shows Rep, the ratio of yields for cen- 
tral and peripheral collisions normalized to the average 
number of nucleon-nucleon collisions in each event sam- 
ple. The ratio is below unity for all pT- The three pt 



regions show different trends as outlined in the discussion 
of Figure Uni (i) In the "soft" region withpr < 2 GcV/c, 
the ratio increases as function oipT- (ii) In the "hard" 
region with px > 4.5 GeV/c, the suppression appears 
to be constant at ~ 0.3, which again indicates that the 
spectra have a similar shape, but with the yield in cen- 
tral collisions being suppressed by a constant factor from 
4.5 to 10 GeV/c. (iii) In the transition region from 2 to 
~4.5 GeV/c, the ratio decreases as function of px- 



B. Suppression of high pr hadrons in Au + Au at 
y/sWN = 200 GeV 

At finite Q^, nuclear modifications of the parton dis- 
tribution (2ft and initial ^i§\ and final state in- 
teractions of the scattering partons can modify the high- 
Pt hadron production rates in hard scattering processes. 
Medium modifications of hadron spectra are often quan- 
tified by the "nuclear modification factor" Raa, which we 
calculate for each centrality class as the ratio of the yield 
per nucleon-nucleon collision in Au + Au to the yield in 
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FIG. 9: (Color online) Ratios of centrality selected pr spectra to the minimum bias spectrum. Ratios for peripheral classes 
are scaled up for clarity. For the pr range shown, most of the systematic errors cancel in the ratio. The remaining systematic 
errors that can change the shape are less than 10% (see Table HVIl and are correlated bin-to- bin in pr- 
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FIG. 10: (Color online) Centrality dependence of (pr 
the average pr of charged particles above a pr threshold as 
defined in Eq. H Shown are {pif""'") values for three p??'" 
cuts, with Pr'" = 0.5, 2 and 4.5 GeV/c respectively. The 
errors shown are statistical only. The systematical errors for 
all data points are less than 3%. 



nucleon-nucleon collisions: 



Raa{pt,ii) = 



1 d^N 



A+A 



coll/ 



N+N 
inel 



dprdrj 



(10) 

{Ncou)/cr^'^i^ is the average Glauber nuclear overlap 
function, (TauAu), for each centrality class. In order to 
calculate Raa , we need a reference spectrum for nucleon- 
nucleon collisions. Due to the lack of charged hadron 
data with sufficient reach in px from our own experiment, 
we construct the N + N reference for charged hadrons 
from the 7r° spectra in p + p collisions at -^/i = 200 GeV/c 
measured by PHENIX [l^, and the charged hadron to 
pion ratio observed in other experiments, as described 



FIG. 11: (Color online) Ratio of charged hadron yields per 
nucleon-nucleon collision between central (0-10%) and periph- 
eral (60-92%) Au + Au collisions. The solid error bars on each 
data point are statistical. The error bar on the left hand side 
of the figure is the overall scale error relative to 0.5, which 
is the quadrature sum of (i) the uncertainty of (Ncoii) (see 
Table and (ii) the uncertainty on the occupancy correction 
(Soccupancy)- The shadcd error band on each data point is 
the Pr dependent systematic error from Sr^i^r^^^^^ and cen- 
trality dependent feed down correction (Sfeeddown) as given 
in Table 1171 



below. 

The PHENIX 7r° spectra from p + p collisions are mea- 
sured out to 14 GeV/c. These data can be parameterized 
by a power-law function. 



1 



" "N+N 



2npT dpxdr] 



= A 



Po 



Po 



(11) 



with A = 386 mb/(GeV/c)2, po 
= 9.99 M- 



Pt , 

1.219 GeV/c, and n 
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In p + p experiments at the ISR, the h/n ratio was 
measured to be 1.6 ±0.16, independent of from 1.5 to 
5 GeV/c, and independent of ^/s from 23 to 63 GeV [3]| . 
Below 1.5 GeV/c, h/ir decreases towards lower pT- The 
ISR data are consistent with data on tt, K, P production 
from FNAL E735 experiment ^=l.S TeV. The 

/i/tt ratio computed from these data increases with pT 
and reaches a value of 1.6 at the end of the measured pT 
range, ~1.5 GeV/c. At high momentum, a h/n ratio of 
~ 1.6 is also observed for quark and gluon jet fragmen- 
tation in e+e~ collisions at LEP by the DELPHI Col- 
laboration [s^- Finally, charged hadron data measured 
by PHENIX \n p + p collisions and data measured by 
UAl vnp+p collisions, both at -^s = 200 GeV/c, give 
consistent h/n ratios when compared to the PHENIX 
p + p tt" data. 

Based on these findings, we assume that /i/tt is con- 
stant above 1.5 GeV/c u\ p -\- p collisions at RHIC and 
that we can scale up the 7r° cross-section (Eq. by this 
factor to obtain a reference for charged hadron produc- 
tion. To be consistent with the data described above, we 
correct this reference below 1.5 GeV/c using an empirical 
function, 



r{pT) 



_ / Rh/7, - a[pn 



\2 
PT) 



for PT < Pmax 
for PT > Pmax 

(12) 

Pmax = 1.6 GeV/c and a — 0.28 
(GeV/c)~^. The charged hadron reference used in this 
analysis is then given by the product of the power-law 
function from Eq. 1111 and the empirical correction from 
Eg.nHas: 



where Rii/tt = 1-6 



1 



^2 h++h- 



2'kpt dpTdr] 



= A 



Po 



Po +PT 



X r{pT) 



(13) 



The systematic errors on the charged hadron N + N 
reference are summarized in Table The main sources 
of uncertainties include: (i) the systematic errors on the 

absolute normalization of the PHENIX tt'' data (<5^orm)i 
which are independent of pT, (ii) uncertainties due to the 

power-law fit to the tt" data {Sj^^f. ) , and (iii) uncertainties 
on Rh/TT (^/i/tt), which are estimated from the spread of 
Rfi/T^ obtained from different data sets used to constrain 

h/TT°. 

Figure El shows the nuclear modification factor 
RAAiPr) for charged hadrons from minimum bias and 
nine centrality classes. The systematic errors on Raa are 
described in the figure captions. At low pT, the charged 
hadron Raa increase monotonically up to 2 GeV/c for 
all centrality classes. At pT > 2 GeV/c, Raa remains 
constant and close to unity for the most peripheral cen- 
trality class. However, in central collisions, it decreases 
at higher pT, down to an approximately constant value 
of 0.2-0.3 for pt > 4-5 GeV/c. This is consistent with 
Figure^2 where the central to peripheral ratio also satu- 
rates above 4-5 GeV/c. This approximately indepen- 
dent suppression pattern has been interpreted as a result 



of the detailed interplay between the Cronin effect, nu- 
clear shadowing, and partonic energy loss [5ll |. 

Also shown in Figure 1121 are Raa for neutral pions 
from ref. The neutral pion Raa values also seem 

to reach maximum around 2 GeV/c, but the changes are 
smaller than those for charged hadrons. Except for the 
most peripheral bin, the neutral pion Raa are always 
below the charged Raa in the range of 2 < < 4.5 
GcV/c. However, at pt > 4.5 GeV/c, Raa for both 
neutral pions and hadrons saturate at roughly the same 
level, indicating a similar suppression for neutral pions 
and charged hadrons at high pt ■ 

The fact that the neutral pion Raa values are smaller 
than inclusive charged hadron Raa at intermediate pt 
{2 < pt < 4.5 GeV/c) has already been observed at 
= 130 GeV [171 . This difference can be explained 
by the large p/ir ratio observed in the same pt range in 
central Au + Au collisions [28L 12^ . This large relative 
proton and anti-proton yield indicates a deviation from 
the standard picture of hadron production at p-r > 2 
GeV/c, which assumes that the hadrons are created by 
the fragmentation of energetic partons. Such a deviation 
has led to models of quark coalescence 33| or baryon 
junctions [s^ as the possible mechanisms to enhance 
the proton production rate at medium pT- Both mod- 
els predict that baryon enhancement is limited to pt < 5 
GeV/c, beyond which jet fragmentation should eventu- 
ally become the dominant production mechanism for all 
particle species. In that case, one would expect a similar 
suppression factor for charged hadron and tt*^, in agree- 
ment with the data at pt > 4.5 GeV/c. Recently, the dif- 
ference of Raa between charged hadrons and pions was 
also argued as the consequence of centrality and particle 
species dependent (fcr) broadening effect |53j . 

If hard-scattering dominates charged hadron produc- 
tion at Pt > 4.5 GeV/c, the particle composition should 
be determined by the jet fragmentation function, similar 
to nucleon-nucleon collisions. Figure [T^ shows h/ir'^ for 
all centrality classes. The systematic errors are explained 
in the figure captions. In the most peripheral collisions, 
the /i/tt" ratio is consistent with the p + p values down 
to ~ 2 GeV/c. In central collisions, the /i/tt^ ratio is 
enhanced by as much as 50% above the p + p value in the 
region 1 < pT < 4.5 GeV/c. This enhancement gradu- 
ally decreases towards more peripheral collisions and re- 
flects the difference of Raa between the charged hadrons 
and TT^s, which is due to large baryon contribution. The 
enhancement also strongly depends on pT'. It reaches a 
maximum between 2.5 and 3.5 GeV/c, then decreases. 
At Pt > 4.5 GeV/c, the /i/tt" ratios for all centralities 
reach an approximately constant value of 1.6, which is 
consistent with the h/n value observed in p + p 'si'l col- 
lisions and in jet fragmentation in e^e~ 32] collisions. 
The similarity of the spectral shape and of the particle 
composition between Au + Au and p + p collisions sug- 
gest that fragmentation of hard-scattered partons is the 
dominant mechanism of particle production in Au + Au 
collisions above pt of 4-5 GeV/c, regardless of the fact 
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TABLE V: Systematic errors on the charged hadron N + N reference spectrum. AU errors are given in percent and are quoted 
as la errors. Positive and negative errors are given separately where appropriate. Most of the errors are correlated with px- 



PT (GeV/c) S]„{%) Sr^,„{%) total (%) 



0.75 
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10.4 


-3.9 


-1- 


9.1 


-15.1 + 5.9 


-18.7 


+ 15.0 


1.00 
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10.4 


-4.1 


-1- 


8.9 


-14.4 + 5.9 


-18.3 


+ 14.9 


1.50 
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10.4 


-4.6 
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8.3 


-11.6 + 5.9 
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+ 14.6 
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-5.1 
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7.7 
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6.7 


-5.9 + 5.9 


-13.3 
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-5.9 + 5.9 


-14.9 


+ 14.3 
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10.4 
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10.5 


-5.9 + 5.9 


-16.0 


+ 15.9 


7.50 
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10.4 


-12.9 


-1- 


14.3 


-5.9 + 5.9 


-17.6 


+ 18.7 


8.50 
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10.4 


-15.8 


+ 


19.4 


-5.9 + 5.9 


-19.8 


+ 22.8 


9.50 


± 


10.4 


-19.3 


+ 


25.9 


-5.9 + 5.9 


-22.7 


+ 28.5 



that the yields do not scale with Ncoii- 

Since Raa values for charged hadrons and tt^s are ap- 
proximately constant at pT > 4.5 GeV/c, we can quantify 
the centrality dependence of the Raa value by calculat- 
ing it from yields integrated above 4.5 GeV/c. The up- 
per panel of Figure [T^ shows Raa for pr > 4.5 GeV/c as 
function of Npart- The Raa values for charged hadrons 
and 7r° agree for all centrality classes within errors. In 
peripheral collisions with Npart < 50, Raa is consistent 
with binary collision scaling. With increasing Npart , Raa 
decreases monotonically, reaching a value of 0.23 ± 0.03 
(0-5% most central) for charged hadrons and 0.24 ±0.02 
(0-10% most central) for 7r°s. There is an additional 14% 
error common to charged hadrons and ir'^s, which origi- 
nates from the uncertainty on the N + N reference and 

NcoU- 

To address suggestions that the yield of high pT 
hadrons in Au + Au collisions may be proportional to 
Npart instead of Ncoii [H, have investigated a 

different ratio, 

Rf-' = 2{N,oa)/ (Npart) X Raa ■ (14) 

^aa"^' for pt > 4.5 GeV/c is shown in the lower panel of 
Figure im together with solid (or dashed) bands repre- 
senting the allowed range if the data follow binary colli- 
sion (or participant) scaling. As discussed above, for pe- 
ripheral collisions, Raa"^* follows more closely the binary 
collision scaling. Above 50 participants, i?^^'"'' varies by 
only ±20%. However, it peaks at (Npart) = 100 and de- 
creases monotonically towards more central collisions 



The decrease of Raa"^* could be a natural conse- 
quence of energy loss of hard scattered partons in the 
medium [H^l- If the energy loss is large, hard scattered 
partons may only escape near the surface of the reaction 
volume. In a cylindrical collision geometry, for which the 
number of collisions from the surface is proportional to 
Npart, binary collision scaling is reduced to an approxi- 
mate participant scaling. Detailed calculations show that 
in this case, i?^^'"'' slightly decreases with Npart depend- 
ing in details on how the energy loss is modelled [s^ . 
This interpretation is also consistent with our previous 
conclusion that, above 4.5 GeV/c, hadron production is 
dominated by hard-scattering although the yield does not 
scale with the number of binary collisions. Gluon satura- 
tion scenarios "2^ also suggest approximate participant 
scaling, with a 30% increase in Raa over the pt range 
4.5-9 GeV/c in central collisions. This increase can not 
be excluded by the data. 



C. Energy dependence and xt scaling 

The inclusive charged hadron and pT spectra and 
h/iT^ ratios suggest that fragmentation of hard scattered 
partons is the dominant production mechanism of high 
Pt hadrons not only in p + p but also in Au + Au colli- 
sions. For p + p collisions this fact was demonstrated on 
general principles well before the advent of QCD by the 
method of "xT-scaling" . This method does not depend 
on whether the initial projectiles are protons or Au ions, 
so it should be directly applicable to Au + Au collisions. 
Since our data show a suppression of high-pr particles 
in central Au + Au collisions with respect to point-like 



^ In the range from 3—4 GeV/c, fi^^"'' for charged hadrons 
is approximately constant, whic h is consistent with earlier mea - 
surements at = 130 GeV |l3 and ^/s]^ = 200 GeV |20| . 

To interpret this constancy as participant scaling is misleading, 
since pion and proton yields change differently with centrality 



in this Pt region, and Rj^'^^* accidentally appears constant for 
inclusive charged hadron. The data above 4.5 GeV/c shown in 
Figure [T^ are free of this effect. 
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FIG. 12: (Color online) Raa for (h'^ + h~)/2 and tt*^ as function of pr for minimum bias and 9 centrality classes according to 
the "Fine" type of centrality classes defined in Table|I] The error bars on the 7r° data points include statistical and systematical 
errors on the Au + Au data and the N + N reference. The error bars on (/i"*" + h~)/2 data points are statistical errors only. The 

common normalization errors {S^orm from Table^ on the references for charged hadrons and tt^s are added in quadrature with 
the uncertainty on (Ncou) and are indicated by the black bar on the left side of each panel. This error ranges from 15% to 36% 
from central to peripheral collisions and can shift all points in the charged and neutral pion Raa up and down together. The 
shaded band on charged Raa includes the remaining systematic errors on the charged N + N reference summed in quadrature 
with the systematic errors from the Au + Au spectra. This error amounts to -12.5% + 18% at low pr and changes to ±12.5% 
at pT = 4.5 GeV/c and ±18.5% at pr = 8 GeV/c. 



16 



scaling from p + p and peripheral Au + Au collisions, it 
is important to investigate whether the production dy- 
namics of high-pT particles in central (and peripheral) 
Au + Au collisions are the same or different from those in 
p + p collisions. We first review the a^^-scaling method in 
p + p collisions and then apply it to the present Au + Au 
data. 

The idea of hard-scattering in -I- collisions dates 
from the first indication of point-like structure inside the 
proton, in 1968, found in deep inelastic electron-proton 
scattering i.e. scattering with large values of 4- 

momentum transfer squared, Q^, and energy loss, v. The 
discovery that the Deep Inelastic Scattering (DIS) struc- 
ture function 



"scales" , or in other words, depends on the ratio 
^ 2Mv 



(15) 



(16) 



independent of as suggested by Bjorken [5g , led to the 
concept of a proton being composed of point-like "par- 
tons" . Since the partons of DIS are charged, and hence 
must scatter electromagnetically from each other in p -I- p 
collisions, a general formula for the cross section of the 
single-particle inclusive reaction 



p + p C + X 



(17) 



was derived |57l | using the principle of factorization of the 
reaction into parton distribution functions for the pro- 
tons, fragmentation functions to particle C for the scat- 
tered partons and a short-distance parton-parton hard 
scattering cross section. 

The invariant cross section for the single-particle in- 
clusive reaction (Eq. I17|l , where particle C has transverse 
momentum pT near mid-rapidity, was given by the gen- 
eral scaling form (53 |: 



E- 



1 



2pj 



where xt = ^pr/^s 



dp^ p'^ \ vs 

(18) 

The cross section has 2 factors: a function F which de- 
pends only on the ratio of momenta, and a dimensioned 
factor, p^", where n depends on the quantum exchanged 
in the hard-scattering. For QED or Vector Gluon ex- 
change [s^l, n = 4. For the case of quark- meson scatter- 
ing by the exchange of a quark [s^, n=S. The discovery 
of high pt pious in p -f p scattering at the CERN-ISR, in 
1972 j5£, ,60, _61 J , at a rate much larger than predicted by 
electromagnetic scattering, but with the scaling form of 
Eq. ^1 proved that the partons of DIS strongly interact 
with each other. 

Inclusion of QCD 62] into the scaling form led to the 
XT-scaling law fEa. [TH|) . 



E 



fa 
dp^ 



1 



G{xt) 



(19) 



where the "a;T-scaling power" n{xT:y/s) should equal 4 
in lowest order (LO) calculations, analogous to the l/q^ 
form of Rutherford Scattering in QED. The structure 
and fragmentation functions, which scale as the ratios 
of momenta are all in the G{xt) term. Due to higher 
order effects such as the running of the coupling constant, 
Q^s(Q^), the evolution of the structure and fragmentation 
functions, and the initial state /cy, measured values of 
n(xT, ^/s) in p + p collisions are in the range from 5 to 
8. 

The compilation of single particle inclusive transverse 
momentum spectra at mid-rapidity from p + p and p + p 
collisions at cm. energy ^/s from 23 to 1800 GeV [sj, 
'so, 63, 64] is shown in Fi gure [T3l for ih^ -I- h~)/2, and 
in Figure 116b for tt" [l^ l65l l66l l67i 16^ ] . The spectra 
exhibit a characteristic shape: an exponential part at 
low Pt < 1 GeV/c which depends very little on ^/s (soft 
physics), and a power-law tail for pt > 2 GeV/c which 
depends very strongly on ^/s (hard physics). The high 
Pt part of the spectra shows a characteristic scaling be- 
havior indicative of fragmentation of jets produced by 
hard-scattering of the quark and gluon constituents of 
the proton as described by QCD [63, IzS 1^ ■ 

The XT-scaling of the single particle inclusive data is 
nicely illustrated by a plot of 



(20) 



as a function of xt, with n(xT,-\/s) — 6.3. The 
(/i+-f /i~)/2 data fFigurelTBb) show an asymptotic power 
law with increasing xt- Data at a given -^s fall below 
the asymptote at successively lower values of xt with in- 
creasing corresponding to the transition region from 
hard to soft physics in the pT range of 1-2 GeV/c. The 
7r° data fFigurelTEb') show a similar XT-scaling but with- 
out the deviation at low xt, since all available data are 
for Pt larger than 1-2 GeV/c. For larger xt > 0.3, a 
value of n = 5.1 |6^ [t^ improves the scaling for the 3 
lower cm. energies, = 38.7, 52.7 and 62.4 GeV. It 
will be a challenge at RHIC to obtain data in this xt 
range to see whether the value of n ~ 5 is the asymptotic 
limit for inclusive single particle production or whether 
n reaches the (LO) QCD value of 4. XT-scaling has also 
been studied in jet production at ^/i = 630 and 1800 
GeV [73, where n = 4.45 is observed in the jet xt range 
0.15-0.3. 

In Au + Au collisions, XT-scaling should work just as 
well as in p + p collisions and should yield the same value 
of n(xT, if high Pt particles are the result of hard- 
scattering according to QCD. This is because the struc- 
ture and fragmentation functions in Au -\- Au collisions 
should scale, so that Eq.^]applies, albeit with a different 
G{xt)- Thus, if the suppression of high-pT particles with 
respect to point-like scaling from p -|- p collisions is due to 
shadowing of the structure functions |43 or gluon satura- 
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FIG. 13: (Color online) Charged hadron to 7r° ratios for minimum bias events and 9 centrality classes according to the 
"Fine" type of centrality classes defined in Table |I] The error bars represent the quadratic sum of statistical and point-by- 
point systematic errors from {h'^ + h~)/2 and tt". The shaded band shows the percent normalization error (dominantly from 
{h+ + h-)/2 data) common to all centrality classes. The dashed line at 1.6 is the /i/tt ratio measured in p -|- p [31| and e'^e 
collisions. 



18 



u 
> 

O 
in 

A 



A ns- 



Au+Au200GeV ■ (h* + h")/2 



> 

0) 

C3 

IT) 
>* 

A 



z < 




FIG. 14: (Color online) Au + Au yield integrated for pr > 
4.5 GeV/c over the N + N yield, normalized using either 

NcoU {Raa in the top panel) or Npart {Raa""* bottom 
panel), plotted as function of (Npart)- The bands represent 
the expectation of binary collisions (solid) and participant 
pair (dashed) scaling. The width of the bands gives the sys- 
tematic errors on Ncou (Npart) added in quadrature with the 
common normalization errors on the N + N references for 
charged hadrons and neutral pions. For charged hadrons, the 
statistical errors are given by the bars. The systematic errors, 
which are not common with the errors for neutral pions and 
which are correlated in pr are shown as brackets. The shaded 
bars around each neutral pion points represent the systematic 
and statistical errors, these errors are not correlated with the 
errors shown for the charged hadron data. 



tion |2g, which are basically scaling effects rather than 
due to a final state interaction with the dense medium, 
which may not scale, the cross sections (Eq. ITi^ at a 
given XT (and centrality) should all exhibit the same sup- 
pression. The initial state shadowing may cause G{xt) 
to change with centrality, but n{xT,^/s) should remain 
constant. In the case of the interaction with the dense 
medium, ZT-scaling may or may not hold, depending on 
the details of the energy loss, for instance, whether or 
not the energy loss of the hard-scattered parton scales 



* There is a slight non scaling effect of the structure functions l47l 
since for fixed xt, changes by a factor of 2.4 between the two 



with its energy. It is also conceivable that the high pT 
particles observed in Au + Au collisions at RHIC have 
nothing to do with QCD hard-scattering [2^ Is^ IH^. In 
this case, striking differences from Eq. E| and the sys- 
tematics observed in p + p collisions should be expected. 

To test XT-scaling in Au + Au collisions, we plot the 
quantities defined by Eq. [201 in Figure El for charged 
hadron and tt" data from ^/snn = 130 GeV and 200 GeV 
for central (0-10%) and peripheral (60-80%) collisions. 
For the power n, we use the same value n(^XT, \/s) — 6.3 
that was used for the p + p data shown in Figure lTSb and 
Figure lTCb . The data are consistent with XT-scaling over 
the range 0.03 < xt < 0.06 for tt^ and 0.04 < xt < 0.075 
for {h+ + h-)/2. 

According to Eq. the ratio of inclusive cross sec- 
tions at fixed Xt equals (200/130)". Thus, the power 
n{xT, \/s) is related directly to the logarithm of the ratio 
of invariant hadron yield at fixed xt as: 



n{xT) 



log{yield{xT, 130GeV)/yield{xT, 200GeV)) 
Zog(200/130) 



(21) 

The power n's for both neutral pions and charged 
hadrons for central and peripheral collisions are shown in 
Figure ITHl While the tt" data in central and peripheral 
collisions and charged hadron data in peripheral collisions 
seem to favor a similar power n, the charged hadron data 
from central collisions require a larger value of n. 

For a more quantitative analysis, the Au + Au data for 
a given centrality and hadron selection are fitted simul- 
taneously for y^SNN = 130 and 200 GeV to the form. 



(xt) 



(22) 



where we have approximated Eg. 1191 bv using a constant 



power n{xT, \/s) and a power-law. 



for G{xt) over 



a limited range in xt- The fit results and errors are 
quoted in Table IVII The corresponding ratios of yields 
are presented by lines in Figure ITHl where the fit ranges 
(0.03 < XT < 0.06 for 7r"s and 0.04 < xt < 0.074 for 
charged hadrons) are indicated by the length of the line. 

For peripheral collisions the fitted values for the power 
are n = 6.33±0.54 and n = 6.12±0.49, for 7r° and charged 
hadrons respectively, which are in quantitative agreement 
with the expectation from p + p collisions. Approximate 
XT-scaling in peripheral Au + Au collisions with the same 
power as observed in p -I- p collisions indicates that hard- 
scattering is the dominating production mechanism for 
high pt particles. In central collisions, neutral pions also 
exhibit XT-scaling with a similar power, n — 6.41 ± 0.55. 
Thus, it seems that high-pT ""^ production is consistent 
with hard-scattering, with scaling structure and fragmen- 
tation functions, for all centralities. 

For charged hadrons, the power found for central col- 
lisions is n = 7.53 ± 0.44. Most of the systematic errors 
are common and cancel between central and peripheral 
collisions, thus the difference of the two powers found for 
charged hadrons. An = Ucent — nperiph = 1-41 ± 0.43 
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FIG. 15: (Color online) (a) Transverse momentum dependence of the invariant cross section at seven center of mass energies 
from different experiments |^|13,E1|H3|. (b) The same data multiplied by ^/s ' , plotted as a function of xt ~ 2pT/\fs- 



compared with that for neutral pion An — 0.09 ± 0.47, 
is significant. 

This difference is consistent with the large proton 
and anti-proton enhancement in central Au + Au colli- 
sions for intermediate pt seen at ^/snn = 130 and 200 
GeV, which appears to violate a^T-scaling. The xt range 
0.04 < XT < 0.074 corresponds to 4 < pT < 7.4 GeV/c at 
^/SNN — 200 GeV, but it corresponds to 2.6 < pr < 4.8 
GeV/c at ^/snn = 130 GeV. If protons are enhanced at 
2 < Pt < 4.5 GeV/c in central collisions at both y^SNN 
= 130 GeV and 200 GeV, then Ucent will be larger than 
nperiph m the measured xt range. Since ^^snn = 200 
GeV data indicate that the proton enhancement is lim- 
ited to the medium pT range, based on the equality of 
Raa for charged hadrons and tt'~' at pt > 4.5 GeV/c 
(Figure J this difference should go away at larger xt ■ 



IV. SUMMARY 



TABLE VI: Results of the simultaneous fit to y'sjvjv ~ 130 
and 200 GeV data using Eq. ^ The fit ranges are 0.03 < 
XT < 0.06 for 71° and 0.04 < xt < 0.074 for charged hadron. 
Only statistical and point-to-point systematic errors on the 
data points are included in the fit, which gives the statistical 
error on n. The normalization errors and other pT correlated 
systematic errors are not included in the fit but are directly 
translated into a systematic error on n. 



Fitting results for over 0.03 < xt < 0.06 


parameters 


0-10% centrality bin 


60-80% centrality bin 


A 


0.973 ± 0.232 


0.843 ±0.3 


m 


8.48 ±0.17 


7.78 ±0.22 


n 


6.41 ± 0.25{stat) 


6.33 ± 0.39{stat) 




±0A9{sys) 


±0.37{sys) 


Fitting 


results for h'^ + h'^ over 


0.04 <XT < 0.074 


A 


2.30 ± 0.44 


0.62 ±0.27 


m 


8.74 ±0.28 


8.40 ± 0.43 


n 


7.53 ±0.18( Stat) 


6.12±0.33(stat) 




±0.40(sys) 


±0.36{sys) 



We have presented a systematic study of the pT and 
centrality dependence of charged hadron production at 
|?7| <0.18 at y/SNN = 200 GeV. The yields per nucleon- 
nucleon collision in central collisions are significantly sup- 
pressed compared to peripheral and nucleon-nucleon col- 
lisions. The suppression is approximately independent 
of Pt above 4.5 GeV/c for all centrality classes, sug- 
gesting a similar spectral shape between Au + Au and 
p + p collisions. At pT > 4.5, charged hadron suppres- 
sion is the same as for neutral pions; the ratio h/ir^ is 



~ 1.6 for all centralities, similar to the /i/tt value mea- 
sured m p + p and e+e^ collisions. The similar spectral 
shape and particle composition at high pT are consistent 
with jet fragmentation as the dominating mechanism of 
particle production in Au + Au collisions for pT > 4-5 
GeV/c. For both charged hadrons and neutral pions, the 
suppression sets in gradually from peripheral to central 
collisions, consistent with the expectation of partonic en- 



20 



1 

cT- -1 
O" 10 
> -2 

S^o 

^ -3 

■g 10 

-4 

■a -5 

D 10 

CO 

■D -6 
UJ 10 

-7 

10 

-8 

10 

-9 

10 

-10 

10 





1 1 1 1 


1 1 1 1 1 








\ (a) 




pbar p 540 GeV: 
UA2 






pp 200 GeV 




. DLJ C M 1 V 

* rnbNIX 






pp 62.4 GeV 






• R110 






CCRS - 






pp 52.7 GeV 


\ ' 




■ R110 J 







□ CCRS 




* 

^ 


pp 38.7 GeV 






' E706 - 








1 




















4- * 


* ^ t t 1 


































t 

. , , , 


1 



10 15 
Pj (GeV/c) 



10 



16 



=^^10 

^10 



10 



11 



10 
10 
10 

io' 

10' 
10 

r, 10 ' 

"> 



CO 

Q. 

■o 

T3 
UJ 



10 
10^ 
10^ 
- 10^ 
10 



C5 



\ 

(b) \ 



pbarp 540 GeV 

» UA2 
pp 200 GeV 

» PHENIX 
pp 62.4 GeV 

• R110 

o CCRS 
pp 52.7 GeV 

■ R110 

□ CCRS 
pp 38.7 GeV 

' E706 



10 



10 




FIG. 16: (Color online) (a) Transverse momentum dependence of the invariant cross section for tt" at five center-of-mass 
energies from different experiments HI El El El El- (b) The same data multiplied by ^/s ' , plotted vs xt = 2pT / \/s- 
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FIG. 17: (Color online) xt scaled spectra for central collisions and peripheral collisions at y^ijvjv = 130 and 200 GeV. The 
left figure shows the tt" xt spectra, and the right figure shows the {h'^ + h~)/2 xt spectra. The central (0-10%) xt spectra 
are represented by triangular symbols, and the peripheral (60-80%) xt spectra are represented by square symbols. The open 
symbols represent xt spectra from ^sjvjv = 130 GeV scaled by a factor of (130/200)®'^. The solid symbols represent xt spectra 
from y/sNN = 200 GeV. The error bars are statistical only. 



ergy loss and surface emission of high px hadrons. xt- 
scaled hadron yields are compared between ^snn — 130 
GeV and ,/snn — 200 GeV Au + Au collisions. We find 
that the xt scaling power n calculated for neutral pions 
in central and peripheral collisions and charged hadron 
in peripheral collisions is 6.3 ± 0.6, similar to p + p colli- 
sions. This again points towards similar production dy- 



namics, i.e. hard-scattering processes as described by 
QCD. However, n is 7.5 ± 0.5 for charged hadrons in 
central collisions, indicating a strong non-scaling modifi- 
cation of particle composition of charged hadron spectra 
from that of p -I- p at intermediate pt, 2-4.5 GeV/c. This 
is consistent with the large /i/tt'^ ratios observed over the 
same pT range in central collisions. 
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FIG. 18: (Color online) The xt scaling power n (according to Eq. I21I I plotted as function of xt calculated for tv" (top-left) 
and [h'^ + h~)/2 (top-right) in central (0-10%) and peripheral (60-80%) collisions. The solid (and dashed) lines indicate a 
constant fit along with the fitting ranges to the central (and peripheral) n{xT) functions. The error bars at each data point 
include statistical and point-to-point systematic errors from y^sjvjv ~ 130 and 200 GeV. The scale errors on xt spectra are 
20.7% (15.9%) for vr'^ xt spectra ratio in central (peripheral) collisions, and 18.6% (15.7%) for (fe^ -I- h~)/2 xt spectra ratio 
in central (peripheral) collisions. These type of errors propagate into the systematic errors on xt scaling power n listed in 
Table ED 
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